THE B-QUADRILATERAL LATTICE, ITS TRANSFORMATIONS 
AND THE ALGEBRO-GEOMETRIC CONSTRUCTION 



ADAM DOLIWA* 



Abstract. The B-quadrilateral lattice (BQL) provides geometric interpretation of Miwa's dis- 
crete BKP equation within the quadrialteral lattice (QL) theory. After discussing the projective- 
geometric properties of the lattice we give the algebro-geometric construction of the BQL ephasiz- 
ing the role of Prym varieties and the corresponding theta functions. We also present the reduction 
of the vectorial fundamental transformation of the QL to the BQL case. 
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1. Introduction 

During last years the many results describing the well known connection between integrable 
partial differential equations and the differential geometry of submanifolds have been transfered to 
the discrete (difference) level, see for example |4"l ll9ll4"3"] . The interest in such a research is stimulated 
from various fields, like computer visualization, combinatorics, lattice models in statistical mechanics 
and quantum field theory, and recent developments in quantum gravity. 
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A succesful general approach towards description of the relation between integrability and geom- 
etry is provided by the theory of multidimensional quadrilateral lattices (QLs) |17| . These are just 
maps x : Z N — > P M (3 < N < M) with planar elementary quadrilaterals. The integrable partial 
difference equation counterpart of the QLs are the discrete Darboux equations (see Section 12.41 for 
details), being found first |Hj as the most general difference system integrable by the d method. It 
should be metioned that the (differential) Darboux equations play an important role in the 
multicomponent Kadomtsev-Petviashvilii (KP) hierarchy, which is commonly considered [121 127| as 
the fudamental system of equations in integrability theory. 




Figure 1 . The geometric integrability scheme 

It turns out that integrability of the discrete Darboux system is encoded in a very simple geometric 
statement (see Fig. 0). 

Lemma 1 (The geometric integrability scheme). Consider points xq, xi, x-i and x% in general 
position in F M , M > 3. On the plane (xQ,Xi,Xj), 1 < i < j ' < 3 choose a point Xij not on the lines 
{xo,Xi), {xo,Xj) and (xi,Xj). Then there exists the unique point x 123 which belongs simultaneously 
to the three planes (x 3} x ls , X23), (x 2 , X12, £23) and (xi, £12, £13)- 

Integrable reductions of the quadrilateral lattice (and thus of the discrete Darboux equations) 
arise from additional constraints which are compatible with geometric integrability scheme (see, for 
example |15II18| ). One of the most important reductions of the KP hierarchy of nonlinear equations 
is the so called BKP hierarchy ^2] (here " B" appears in the context of the classification theory of 
simple Lie algebras) . In [33] it was shown that the r- function of the BKP hiererchy satisfies certain 
bilinear discrete equation (the transformation between the infinite sequence of times of the hierarchy 
and the corresponding discrete variables is called the Miwa transformation) 

(1-1) T r (123) = r (12)T(3) - T (13) r (2) + T(23)T(i), 

which is known as the discrete BKP or the Miwa equation. Here and in all the paper, given a fuction 
F on 1 N , we denote its shift in the ith direction in a standard manner: Fu\{ji\, . . . , rij, . . . , tin) = 
F(m, . ..,rii + 1, .. .,njv). 

The linear problem and the Darboux-type (Moutard) transformations for the discrete BKP equa- 
tion Hl.lfl were constructed in |38| . In literature there are known several geometric interpretations 
of the discrete BKP equation in terms of the reciprocal figures and inversive geometry |29l l4~3"] , or in 
terms of the trapezoidal nets 0j . It should be also mentioned that the discrete BKP equation has 
been recently investigated, under the name of the cube recurrence, in combinatorics |4()l 1251 18] . 

In this paper we propose (see Section |2J) another geometric interpretation of the discrete BKP 
equation, which we consider from the point of view of the quadrilateral lattice theory. This new 
reduction of the quadrilateral lattice, which we call the B- quadrilateral lattice (BQL), is projectively 
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Figure 2. Elementary hexahedron of the B-quadrilateral lattice 



invariant and is based on additional local linear constraint. Sectional is of rather elementary geo- 
metric nature, but the results obtained there have far reaching consequences. In fact, the paper 
gives new arguments supporting the conjecture that basic integrability features are consequences of 
incidence geometry statements (see also introductory remarks in 

More involved techniques are used in Section where we elaborate the algebro-gemetric method 
to produce large classes of the B-quadrilateral lattices and the corresponding solutions of the discrete 
BKP system of equations. In doing that we start from recent results of [221 > where restrictions on the 
algebro-geometric data of the discrete Darboux system £Q compatible with the discussed reduction 
were given. Then we proceed to formulas for the wave and r-functions of the BQL in terms of the 
Prym theta functions related with the algebraic curves used in the construction. We transfer this 
way the algebro-geometric method of construction of solutions of the BKP hierarchy |2j and of its 
two-component generalization 07] to the discrete level. 

We also present, in Section 0] the corresponding reduction of the vectorial fundamental transfor- 
mation of the quadrilateral lattice |2J and establish its link with the Pfaffian form of the vectorial 
Moutard transformation found in |38| . In Appendix we give alternative proof of a crucial auxilliary 
result of the paper, and we summarize basic properties of Pfafhans. 



We start with discussing the gometric constraint, which imposed on the quadrilateral lattice 
allows to define its new integrable reduction. Then we proceed to the algebraic description of such 
reduced lattice showing its connection with the discrete BKP equation. Finally, we discuss relation 
between the r-function of the B-quadrilateral lattice, and the r-function of the quadrilateral lattice. 

2.1. Geometric definition of the BQL. 

Proposition 2. Under hypothesess of Lemma^\ assume that the points xq, x\i, x\z, X23 are 
coplanar, then the points x\, x%, X3, and X123 are coplanar as well (see Figure^). 

It can be shown by the standard linear algebra (for a synthetic- geometry proof see a remark below). 
We perform however the calculations, because the way we are going to do it will be important in 
next Sections in showing connection of the BQL with the discrete BKP equation. 

Lemma 3. Under hypotheses of Proposition^ for fixed initially homogeneous coordinates Xq and 
Xi (gauges) of Xq and x\, there exist a gauge such that the following linear relations hold 



2. The B-quadrilateral lattice 



(2.1) Xij -x = P{ Xl - Xj ), l<i<j<3, 

where the coefficients f lJ depend on the actual positions of the points Xij . 



Proof. The coplanarity of the four points Xo, x\, x 2 and x\2 can be algebraically expressed as the 
linear relation 

ctxo + (3xi + jx 2 + 6x 12 = 0, 
where, by the genericity assumption (no three of the points are collinear), all the coefficients do not 
vanish. By plaing with rescalling the homogeneous coordinates of x 2 and x\ 2 we can transfer above 
equation to the form [j2.1J) 

(2.2) x 12 -x Q = f 12 {x 1 -x 2 ). 

Similarly, we can rescale the homogeneous coordinates of X3 and X13 to express planarity of the 
corresponding elementary quadrilateral as 

(2.3) x 13 -x a = f 13 {x 1 -x 3 ). 

However, with fixed gauges Xo, x 2 and x 3 the coplanarity of Xo, x 2 , x 3 and x 2 3 can be expressed, 
by plaing with the gauge of x 2 %, at most as 

(2.4) x 23 - x = ax 2 - bx 3 . 
Then 

(2.5) x A X12 A x 13 A x 23 = f 12 f 13 (a - b)x A x\ A x 2 A x 3 , 

and at this moment we use coplanarity of xo, x 2 , x\ 3 , x 23 , which is equivalent to a — b. □ 

Remark. Notice that the whole reasoning can be applied even if the gauges of Xo and of X\ are not 
fixed initially (we could then achieve f 12 = 1). However we will need this additional restriction in 
next Sections. 

Proof of Proposition^ By the linear algebra, the homogeneous coordinates of the point x\ 23 £ 
(£3, £13, x 23 ) n {x 2l x 12 , x 23 ) fl {xi,xi2, £13), in the gauge of Lemma |3] read 

= (f 13 f 12 + X! - (V 3 + r - ^) X2 + f/ 13 - z 23 + X3 , 



p~ L ^ \ J J ' J23 / 1 W ■> J13 I ' I ' ' f 1 

(we still keep the undetermined yet factor p) which gives coplanarity of x\ 23l x\, x 2 and x 3 . □ 
Corollary 4. By fixing the gauge function p to 

f 13^23 j-12 j-13 fl2f23 
_ _ 9 j-13 _ 9 j-12 _ --23 , J J 1 J J 1 J J 
P — Z J Z -l J + J12 + f23 + fl3 ' 

we find that the linear relations on the new facets (containing x± 23 ) of the cube are of the form (|2.1|l 
again, for example 

xi 23 -xi = f 23 (x 12 - X13), 

where 

f23 _ I 

Jl J12 jl3 „ J12 p3 _|_ J13 J23 ' 

Remark. For a geometrically oriented Reader we would like to comment on another interpretation 
of Proposition |21 It is related to the notion (see, for example ^Hl) 01 the quadrangular set of points 
which are the intersection points of the lines of a complete quadrilateral (add the diagonals). Such a 
configuration is usually denoted by Q(ABC, DEF), where the first three points A,B,C lie on sides 
through one vertex while the remaining three D, E, F lie on the respectively opposite sides, which 
form a triangle. It is known that Q(ABC, DEF) implies Q(DEF, ABC). In notation of Proposition 
[21 denote by £ the intersection line of the plane (xi, x 2 , x 3 ) with the plane (x\ 2 , x 23 ,x\ 3 ) (containing 
also the point xo). Denote by A, B, C, D, E, F intersections of sides of the complete quadrilateral 
with vertices xq, x\ 2 , x 23 ,xi 3 with £ (see FigureEJ, i.e. Q(DEF, ABC). The statement of the Lemma 
is equivalent to the fact that the lines {A,x\), (B,x 2 ), (C,x 3 ) intersect in one point (which is X123); 
see Excercise 1 of Section 2.4 of fTT)| . 
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Figure 3. Two quadrangles 



Remark. As it was pointed to me by Yuri Suris, Proposition's equivalent to the Mobius theorem 
on mutually inscribed tetrahedra. Indeed, vertices xo, x\, X2, X3 of the tetrahedron {xo, x%, X2, X3} 
are contained in the facial planes of the tetrahedron {2:12, X13, 2:23, 2:123}, and vice versa. In fact, 
Figure appears in Mobius' original proof of the theorem. 

Remark. Proposition can be also considered as a special version of the Miquel theorem jSH] (used 
in [Hj to show integrability of the circular lattice) in the same way like Pappus' hexagon theorem is 
a special case of the Pascal theorem. 

We conclude this Section by defining new reduction of the quadrilateral lattice. 

Definition 1. A quadrilateral lattice x : Z, N — > P M is called the B- quadrilateral lattice if for any 
triple of different indices k the points x, xuj\, xufy and xuja are coplanar. 

Corollary 5. In the B- quadrilateral lattice, for any triple of different indices i,j,k the points Xu-\, 
X(j) , X(fc) and X(ijk) are coplanar. 

2.2. Multidimensional consistency of the BQL constraint. As it was shown in ^7] the pla- 
narity condition, which allows to construct the point £123 as in Lemma ^ does not lead to any 
further restrictions if we increase dimension of the lattice. This is the consequence of the of the 
following geometric observation. 

Lemma 6. Consider points xo, x\, X2, X3 and X4 in general position in ¥ M , M > 4. Choose 
generic points Xij S {xq, Xi, Xj), 1 < i < j < 4, on the corresponding planes, and using the planarity 
condition construct the points Xijk € {xq, Xi, xj, Xk), l<i<j<k<4 - the remaining vertices of 
the four (combinatorial) cubes. Then the intersection point X1234 of the three planes 

(£12,2:123,2:124}, (2:13, £123, 2-134), (2:14,2:124,2:134) in (2:1,2:12,2:13,2:14), 
coincides with the intersection point of the three planes 

(2:12,2:123, 2:124), (2:23,2:123,2:234), (2:24,2:124,2:234), in (2:2, 2:12, 2:23, 2:24), 
which is the same as the intersection point of the three planes 

(2:13,2-123,2:134), (2:23,2:123,2:234), (2-34, 2:134, 2:234), in (2:3,2:13,2:23,2:34), 
and the intersection point of the three planes 



,2:14,2:124,2:134;, 



'2:24,2:124,2:234), (2-34,2:134,2:234 



in 



X4,Xu,X24,X3i). 



Remark. In fact, the point X1234 is the unique intersection point of the four three dimensional 
subspaces (xi, x\ 2 , #13, ^u), (%2, %12, 2:23, ^24), {#3, £13, 2:23, £34), and (x 4 , ecu, ^24, 2:34) of the four 
dimensional subspace (xo, Xx, X2, 2:3, X4). This observation generalizes naturally to the case of more 
dimensional hypercube with the planar facets. 

The goal of this Section is to show an analogous result for B-quadrilateral lattice. Notice, that in 
previously known reductions of the quadrilateral lattice, such as the symmetric ^H] or the quadratic 
[T5] reduction, the additional constraint was imposed on initial quadrilaterals. Then the multidi- 
mensional consistency of the reduction was the result of the three dimensional consistency of the 
constraint and the multidimensional consistency of the quadrilateral lattice. 

In the BQL case, the constraint is imposed on the level of elementary cubes. Therefore its four 
dimensional consistency is crucial for integrability of the B-quadrilateral lattice, and once proven, 
implies consistency of the reduction in more dimensions. 

Proposition 7. Under hypotheses of Lemma\Q assume that the BQL condition holds for the initial 
data, i.e., the point x$ belongs to the four planes (xij,Xik,Xjk), l<i<j<fc<4. Then all the 
three dimensional (combinatorial) cubes obtained in the construction satisfy the BQL constaint, i.e., 

X\ £ (£123, £124, Zl34}, ^2 S (£123, a;i24,a:234), 2:3 £ (£123, 2;i34, X234}, 2:4 £ (X124, 2:134, X 2 34)- 

Proof. Consider the gauge of Lemma |3J If we add into the construction points 2:4 and X14,, then by 
fixing suitably their gauges X4 and X14, we can rewrite the coplanarity condition of xq, x%, X4 and 
Xu in the form l|2.1[l . The same argument like in the proof of Lemma |3] implies that the algebraic 
coplanarity conditions of xq, Xi, 24 and x^, i — 2, 3 take the form of equation (|2.1(l . 
By fixing gauges of points x^k like in Corollary^ we obtain the relations 

(2.6) Xij k - ^ = f-(xij - Xik), i,j,k disctinct, 



fijfik _ fij fjk _|_ fikfjk ' 
with P l = -/«. 

In equations H2.6(l let us fix i = 1 and consider the three pairs (J, k): (2, 3), (2, 4) and (3, 4). Then 
after simple calculation we obtain the following linear relation 

(2.8) A 24 A 34 (*i23 - *i) - /i 23 /i 34 (*i24 - xi) + / 1 23 / 1 24 (a; 1 34 - xi) = 0, 

which shows that x% € (2:123, 2:124, 2:134). Other cases are similar. □ 

Corollary 8. Under assumptions of Proposition ^ the point X1234 belongs to the four planes: 
(X12, 2:13,2:14), (X12, 2:23, X24), (2:13,2:23,2:34) and (xu, x 2 4, x 3 a)- 

Remark. The same procedure can be applied when we increase dimension of the hypercube keeping 
the BQL constraint. 

2.3. BQL and the discrete BKP equation. 

Proposition 9. A quadrilateral lattice x : *Z N — > ¥ M is a B-quadrilateral lattice if and only if it 
allows for a homogoneous representation x : 7L N — > R* /+1 satisfying the system of discrete Moutard 
equations (the discrete BKP linear problem) 

(2.9) x (ii) -x = f 3 (x (l) -x u) ), l<i<j<N, 
for suitable functions f % i : Z w — > R. 

Proof. As we have shown above (we present an alternative difference-equation theory proof in the 
Appendix), the B-quadrilateral lattices indeed allow for such a gauge; here the Remark after LemmaOl 



where 
(2.7) 



turns out to be important. Conversely, three equations (|2.9[1 for the pairs (i,k), (j,k) imply 

the linear relation 

(2.10) f jk f ik (x m -x) + f ij f ik (x (jk) -x)- f ij f jk (x m -x) = 0, l<i<j<k<N, 
expressing coplanarity of the four points x, Xuj\, xuy\ and Xuf.y □ 

The system (I2.9|) is well known in the literature |2H]. Its compatibility leads to the following set 
of nonlinear equations 

(2.11) 1 + /*(/«_ /*) = /*/« = i,j,k distinct, 
with p l = -fi. 

Remark. The system <|2.10[) can be actually solved, see for example 0], (we have used this fact in a 
hidden form already). Simply replace the lower index i in the "star-triangle" relation 1|2.7() by the 
shift (i). 

On the other side, the second equality in the compatibility condition implies existence of the 
potential t : Z N — > R, in terms of which the functions / y can be written as 

(2.12) f v = mm t 

The first equality can be then rewritten in the form of the system of the discrete BKP equations 

(2.13) TT (i 3 k) = T (i 3 ) T (k) ~ T (ik)T(j) + T(jk)T(i), 1 < i < j < k < N. 

Remark. Two dimensional quadrilateral lattices whose homogeneous coordinates satisfy (up to a 
gauge) equation (|2.9|) are characterized geometrically [23 by condition that any point x and its 
four second-order neighbours £(±i±2) a- re contained in a subspace of dimension three. Obviously, 
any two dimensional slide of the B-quadrilateral lattice fulfills this property, which can therefore 
serve as definition of a two dimensional BQL. However the example of the standard injection Z w — > 
R N C P N shows that without additional requirements this property does not characterize completly 
multidimensional BQL. 

2.4. The r-functions. In this Section we present relation between the r-function of the quadrilat- 
eral lattice ^Hj; which we donote here by f , and the above r-function of the B-quadrilateral lattice. 
From the relation between the KP and BKP hierarchies we expect that within the class of the 
B-quadrilateral lattices f should be equal to the square of r. 

Let us recall briefly the algebraic construction of the T-function of the quadrilateral lattice (the 
geometric meaning is presented in [IB])- The non-homogeneous coordinates x : 7L N — > R M (we 
restrict our attention to the affine geometric aspects of the theory) of the quadrilateral lattice 
satisfy the system of Laplace equations 

(2.14) x {lj) - x = a l3 (x {l) - x) +a rt (x u) - x), i^j. 

The functions a u : Jj N — > R are not arbitrary (the system (|2.14|) must be compatible) , in particular 
they can be parametrized in terms of the potentials h l (the Lame coefficients) as follows 

(2.15) o« = ^, 

Define the so called rotation coefficients /Jy from equations 

(2.16) Aihj = h^pij, i ^ j, 
and the normalized tangent vectors Xi from 



(2.17) 



AiX — hiXi. 



Then the Laplace system (|2.14[) takes the first order form 

(2.18) & j X i =0 fj X j , i^j, 
and its compatibility reads 

(2.19) Ajfiik = f3ij(k)/3jk, i,j,k distinct. 

The discrete Darboux equations (|2.19|l imply existence of the potential f (the r-function of the 
quadrilateral lattice) 

(2.20) Zl^L = l-^ ij0jh i^Lj. 

T (i) T U) 

In looking for the Lame coefficients hi in the reduction from QL to BQL we can compare both 
linear systems (|2.9|l and 1)2.14(1 . and the expresions l|2.15|l and (|2.12|l to obtain 

(2.21) h = (-l)E — . 

The corresponding rotation coefficients are then given by (below we assume i < j) 

(2.22) /Jy = -(-l)Ei<*<i m * 
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(2.23) /3 Ji = -(-l)Si<*<j'»* 
which implies (compare with formula H2.20[l ) 

(2.24) 1 - A,/^ - f , MJ- 
Therefore, we can summarize the above considerations as follows. 

Proposition 10. Given B- quadrilateral lattice x with the r-function r. Then, formally on the level 
of the reduction of the system of discrete affine Laplace equations (|2.14|) to the system of discrete 
Moutard equations (|2.9|l . the Lame functions and the rotation coefficients are given in terms of t by 
equations 12.21|I - H2.23(I . and the corresponding r-function f of the quadrilateral lattice is given as 

(2.25) f = r 2 . 

Remark. We should be aware that the r-function of the quadrilateral lattice is defined with respect 
to the coefficients of the affine Laplace equation. Equation l|2.9|l . although formally written in the 
affine form, is a consequence of the projectively-invariant definition of the B-quadrilateral lattice. 
Therefore, in this formal correspondence the geometric meaning of the rotation coefficients in the 
BQL reduction has been lost. We mention that the affine geometric meaning of equation 1)2. 9[) . and 
therefore also of the rotation coefficients l|2.22|l - (|2.23|l . can be provided within the context of the 
trapezoidal lattices 0] (see also 03 Ell) ■ 



3. Algebro-geometric construction of the BQL 

Below we apply the algebro-geometric approach, well known in the theory if integrable systems 
to the B-quadrilateral lattice reduction. It is known (221 that in the BQL reduction case, the 
generic algebraic curve, used to generate solutions of the discrete Darboux equations pQ, should be 
replaced by a curve admitting a holomorphic involution with two fixed points. Such curves have 
been already used in construction of solutions of the BKP hierarchy |14j and of its two-component 
generalization 07] . We develop the corresponding results of [221 an d we present the explicit formulas 
for the lattice points and the solutions of the discrete BKP equation in terms of the Prym theta 
functions related to such special curves. 



3.1. Curves with involution and their Prym varieties. Let us first summarize some fact from 
theory of Riemann surfaces (see j2U Consider T A T a ramified double covering of genus 

g = 2g of a compact Riemann surface T of genus g with exactly two branch points Qo ; Qoo- Denote 
by a : T — > T the holomorphic involution permuting sheets of the covering, i.e., T = T/a. 

The map it* : J{Y~) — > J(X) lifting divisor classes of degree is an injection. The holomorphic 
involution a extends to J(T) and allows to define the Prym variety 

(3.1) V a (f) = {A-a(A)\AeJ(f)}. 

The natural epimorphism i : J{V) x V(r(T) — > J^(T) has finite kernel consisting of 4 9 half-periods in 

Jin 

There exists a basis of cycles <2j, &i, 1 < i < 2<? on T with the canonical intersection matrix such 
that 7r(afe), it(bk), 1 < A; < <?, is a canonical basis of cycles on T and 

a(a k ) = -a g+kl a{b k ) = -b g+kl 1 < k < g. 

The corresponding normalized holomorphic differentials u>i 

(3.2) I tdi = 5 tj , l<i,j<2g, 
satisfy 

(3.3) <j*(uo k ) = -u g+k . a*(u g+k ) = -w fe , 1 < /c < 3. 
The differentials 

Ufc = uj k - uj g+k , c*(u fe ) =u k , 1 < fc < g, 
form a basis of normalized holomorphic differentials on T, while the odd differentials 

(3.4) w k = Lo k + ujg +k , a*{w k ) = ~Wk, 1 < fc < g, 

are called normalized holomorphic Prym differentials. Then the Riemann matrix 

(3.5) Bjk = f u k , 1 <j,k <g 
for f 1 has the form 

(3 - 6) B= 2{n-B n + B 

where 

(3.7) B lk =<fu k , l<j,k<g 

Jbj 

is the corresponding Riemann matrix for T, and II is the matrix of the 6-periods of the Prym 
differentials 

(3.8) Il jk = I w k , l<j,k<g. 

The matrix II is symmetric and has positively defined imaginary part, and defines the Prym theta 
function 6>(z;IT), z e C 9 , 

(3.9) 6>(z; n) = ^2 ex P {^(i, nn) + 2vri(n, z)} , 

where (•, •) denotes the standard bilinear form in C s . 

With the above choice of the period matrix B for T and with Qoo as the base-point of the Abel 
map 

(3.10) A(P) = f u, u = (u 1 ,...,u 2g ), 



the lift of a from .7(f) = £ 29 /{I 2g , B) to C 2g reads 

a(Zl, . • • , Z 2 g) = — . . . , Z 2 g, 21, ... , Zg), 

while the map 7r* : <J(T) — -> J7(r) is represented by 

7T*(Zl, • ■ ■ ,Zg) = {Zl,...,Zg,-Zl,...,-Zg). 

The Prym variety is a principally polarized abelian variety isomorphic to V = C 9 / (Ig,H), and the 
injection : C 9 /(I g , II) -> C 29 /(I 2g ,B) is given by 

0(21, . . . = (zi,...,Z ff ,Zl,...,Z g ). 

There holds the following analog of the Riemann theorem. 

Lemma 11 (|2SI)- Given e G C ff suc/i iftat (e; II) 7^ 0, i/ien i/ie zero divisor Z of 9 ^ J^" 10 — e; II 

on r is of degree 2g and satisfies the relation 

(3.11) <Xe) = A(Z)- A(Q )+7T*K r , mod(I 29 ,fl), 

where Kr is i/ie Riemann constants vector of V . Moreover Z + a(Z) — Qq — Qoo is a canonical 
divisor on T. 

By ujst, (S, T G r, S T) denote the unique meromorphic differential holomorphic in T \ {S, T} 
with poles of the first order in 5, T with residues, correspondingly, 1 and —1, and normalized by 
conditions 

(3.12) <f uj st = 0, l<j<2g. 

J a,j 

It is known that the o-periods of such differentials are given by 

(3.13) f uj ST = 2Tri Uj, l<j<2g, 

Jbj JT 

with the integral being taken along a curve joining S to T in f \ \J 2 / =1 Oj \ Uj£i bj. Moreover, 
the following relation between two such differentials holds (with the paths of integration being 
appropriately choosen [21]) 

(3.14) J ujst = J u P q. 

Remark. Such a form can be expressed in terms of the theta function on T as 

0(A(P)-A(S)-frB) 

WST = CtP log 5 — , 

where £ is a general point of the divisor @ of zeros of the theta function in J(T). 

3.2. Explicit algebro-geometric formulas. Given N points Qi G T, different from Qo,Qoo, and 
an effective non-special divisor D of degree 5 such that 

(3.15) D + a{D)-Q -Q oo J2i, = (r) C t 

is a canonical divisor. For an arbitrary m G Z w there exists |22| the unique function ip{m) mero- 
morphic on r having in points Qi (in points o~{Qi)) poles (corrspondingly, zeros) of the order m„ 
no other singularities except for possible simple poles in points of the divisor and normalized to 
1 at Qoo. In |22) it was shown that, as a function of the discrete parameter m, the wave function ip 
satisies the system of the discrete Moutard equations 

(3.16) ij (ij) (P)-TP(P) = fi(iP (i} (P)-i> {j )(P)), l<i<j<N, 



where 

(3.17) r = Hm = lim ^ ^ , i < J. 

To obtain B-quadrilateral lattices we pick up M + 1 points Pi, . . . ,Pm+i G T of the Riemann 
surface. Then Xj(m) — \&(m|P,), 1 < j < M + 1 serve as homogeneous coordinates of the lattice. 
However in this way we obtain B-quadrilateral lattices in complex projective space. In order to get 
real lattices certain additional restrictions, which were also given in [22], should be imposed on the 
algebro-geometric data. 

In |22] the multidimensional aspects of the system (|3.16|l were not of particular importance. Also 
the role of the Prym variety and of the corresponding theta function was not fully exploited. Our 
goal here is to fill this point. We start from the imediate consequence of Definition (|3.1|l of the Prym 
variety. 

Corollary 12. Denote by D{m) the divisor of additional zeros ofi/)(m), then 

. N 

(3.18) D(m) - D V m k (a{Q k ) - Q k ) 



fc=i 



moves linearly within the Prym variety. 



An important part of the algebro-geometric theory of integrable systems consists on providing 
the explicit formulas, in terms of the Riemann theta functions of the corresponding Jacobi varieties, 
for the wave functions and the soliton fields. In the case of the special Riemann surfaces used in the 
paper, there exist [23 formulas connecting the theta functions of T, T and V a - However, instead of 
reducing the explicit expressions given in [Tl I22j for the generic curves, we will follow the reasoning 
of |14) . In order to present the explicit formulas, in terms of the (Riemann-) Prym theta function, 
for the wave function and other relevant data we will use Lema ITT1 

Let us define 

rQk 

V fc = / weC 9 , 1 < k < N, w = (w 1 ,...,w g ), 



then equations (j3.4|l and H3.13(l imply that 

(3.i9) w^r^-^f^, b=(h,...,b 2g ). 




<r(Q k ) 2ni Jb 

Proposition 13. The BQL wave function tp(m) can be written down with the help of the Prym 
theta functions as follows 

(3.20) ip(m P) = ^— = exp 

*(E£Li m k V k + e; H) 8(f^w - e; II) 

where 

(3.21) 0(e)=A(-D)-A(Q o ) + 7r*(Kr). 
Proof. Using the property (|3.15() of the divisor D, the Hurwitz formula 

(3.22) C t j2 ^ {T) n*C T + Qo + Qoo 
relating canonical divisors on T and T, and the relation 

(3.23) A(C r ) = -2K r 
between the canonical divisor and the Riemann constants vector, we obtain 

a(A(D) - A(Q ) + 7r*(K r )) = -A(D) + A(Q„) - vr*(K r ), 
which asserts that the definition of the vector e in Ij3.21|l is meaningful. 



To show that the right hand side of equation i|3.2U[l is single valued on f we check that it is inde- 
pendent on the integration path in the integrals f„ . When two paths differ by an elementary cycle 
we use the properties Ij3.2|l - H3.8jl 01 the holomorphic differentials, the quasi-periodicity properties of 
the theta functions 

0(z + e fe ; n) = 6(z; II), 9(z + He fe ; n) = cxp(-^n fcfc - 2mz k )6{z] II), 

where efc are vectors of the standard basis in C 5 , and the relations Ij3.12|l and (|3.19|1 . From now on 
our path of integration avoids the cuts, like in fromulas 3|) - (I3.14|) . 

As the normalizaton condition at Qoo is obvious (the theta function is even) we are left with the 
analyticity properties. Lemma ITTI implies that the right hand side has simple poles at points of the 
divisor D. Apart from the zeros of the theta function in the nominator (which may eventually cancel 
with the poles at D), the only other poles and zeros are consequences of the analytical properties of 
the integral in the exponential part. Let us choose a local parameter z k (P) at Q k , then 



which implies that 

r ^<r(Q k )Q k ^ -l0g* fc (P) + 0(l), 

and, in consequence, the right hand side in equation (|3.20|) has pole of order m k at Q k . Similarly, 
since z k (a(P)) is a local parameter at cr(Q k ), we have 

"*{Qu)Q* P ^= Qk) logz k (a(P))+0(l), 
and the right hand side in equation (|3.2()(1 has zero of order m k at a(Q k ). O 



Corollary 14. The potentials read 

g(v« + ££L 

m kV k + e; fl) 9(\i + V, + ™feV fc + e; II 



(3.24) f J {m) = ■ 77 ^ N — — ■ —Ay , < J 



where 

(3-25) Xij = exp (j^ uj aiQj}Q ^j , i < j. 

The BQL (the discrete BKP) T-junction within the above class of solutions reads 

N 

(3.26) r(m) = 0(£ m k V k + e; IT) J] A™ . 

k =l i<j 

Proof. Expression Ij3.24|l for the potentials is a direct consequence of their algebro-geometric defi- 
nition l|3.17fl and of equation l|3.20fl . Then the formula l|3.26f) for r-function follows easily from its 
defmiton lj!TT2j) . □ 

Remark. From general considerations of |22j we know that 

(3-27) = exp ( / uj (7{Q: . ) q j ) = - exp ( f ^(q^qA , i < j, 

\JQao ) \JQ^ ) 

which reflects the second equality in (|3.17(l . 



4. Transformations of the B-quadrilateral lattice 



Below we present the reduction of the vectorial fundamental transformation compatible with the 
B-quadrilateral lattice constraint. In literature there is known the direct vectorial Moutard 
transformation between solutions of the BQL linear problem 1)2. 9)1 providing thus the corresponding 
transformation between solutions of the discrete BKP equation 1)2.13(1 . Our goal will be to find the 
transition to the Pfaffian expressions of starting from the BQL reduction of the fundamental 
transformation. In describing this connection we follow the ideas of where similar problem 
between the Grammian expressions for binary Darboux transformation of the KP hierarchy has 
been transformed, in the BKP reduction, into the Pfaffian form [2H| (see also QUI for other 
aspects of the relation of Pfaffians with the BKP hierarchy and the discrete BKP equation). 

4.1. The fundamental transformation of the QL. Let us first recall some basic facts concerning 
the vectorial fundamental transformation of the quadrilateral lattice. Geometrically, the (scalar) 
fundamental transformation is the relation between two quadrilateral lattices x and x such that for 
each direction i the points x, x, xu) and xu) are coplanar. 

We present below the algebraic description of its vectorial extension (see jHSIISESI for details) in 
the affine formalism. Given the solution Y i : Z N — > V, V being a linear space, of the linear system 
(j2~T£|) . and given the solution Y* :Z N ^V*,V* being the dual of V, of the linear system fFHty . 
These allow to construct the linear operator valued potential tl(Y, Y*) : Z N — > L(V), defined by 

(4.1) A i fl(y i Y') = Y i ®Yl i=l,...,N; 
similarly, one defines Sl(X, Y*) : Z N -> L(Y, R M ) and Sl(Y, h) : Z N -> V by 

(4.2) A i Sl(X,Y*)=X i ®Y*, 

(4.3) Aifl(Y,h) = Y t ® hi. 

If Sl(Y, Y*) is invertible then the (vectorial) fundamental transform of the lattice x is given by 

(4.4) x = x- Sl(X, Y*)fl(Y, Y*y 1 n(Y, h). 

The corresponding transformation of the r-function f of the quadrilateral lattice reads 

(4.5) f = fdetn(r,y*). 

When dimV = 1 we obtain the formula relating the quadrilateral lattice x and its fundmental 
transform x. The vectorial fundamental transformation can be considered as superposition of dim V 
(scalar) fundamental transformations; on intermediate stages the rest of the transformation data 
should be suitably transformed as well. Such a description contains already the principle of per- 
mutability of such transformations, which follows from the following observation |21| . 

Lemma 15. Assume the following splitting of the data of the vectorial fundamental transformation 
(4-6) Y i =(*{), Y* = (Y* ai , Y* H ), 

associated with the partition V = V a © Vt>, which implies the following splitting of the potentials 
(4 7) MY h) ( n{Ya > h) \ n<Y Y*) ( n{Ya ' Y:) n{Ya > Y * b) \ 

(4.8) n(x, y*) = ( n(x, Y* a ) , n(x, Y* b ) ) . 

Then the vectorial fundamental tansformation is equivalent to the following superposition of vectorial 
fundamental transformations: 

1) Transformation x -> x {a} with the potentials fl(Y a : h), ft(Y a ,Y* a ), fl(X,Y* a ) 

(4.9) £ {a} = x - n(X, Y*MY a , Y*)~ 1 £l(Y a , h). 



2) Application on the result the vectorial fundamental transformation with the transformed potentials 



(4.10) x {a ' b} = x {a} -n(X,Yt)^[Cl(Y b ,Yl)^}- 1 Cl(Y b ,h) {a} , 
where 

(4.11) Ci(Y b , h)W = n(Y b , h) - n(Y b , Y* a )n(Y a ,Y* a )- 1 n(Y a , h), 

(4.12) tl(Y b , Y* b )W = n{Y b , Y* b ) - n(Y b , Y* a )n{Y a , YD^niY 11 , Yl), 

(4.13) n(x, n) {a} = n(jc, n) - n(jr, r;)n(y a , nr 1 ^*" 1 , n)- 

Corollary 16. TTie normalized tangent vectors Xi and the Lame coefficients hi are transformed, 
at the intermediate step, according to formulas 

(4.14) x\ a} = x t - n(x, Y* a )n(Y a , Y* a )- l Yi, 

(4.15) h\ a] = hi- Y* a n(Y a , Yiy'niY", h), 

which also give the corresponding transforms of the second set of transformation data Y b and Y b 

(4.16) Y f a} - Y b - n{Y b , Y* a )n(Y a , Y* a )~ l Y a iy 

(4.17) Y tb ia} = Y* b Y* a Cl(Y a , Y* a )~ 1 Sl(Y a , Y b )] 
which agree with the transformation rules (14.13(1 for the potentials, i.e., 

n(Y b ,h)^ = n(Y {a \h^), 



b h 

n(x,Yi)^ = n(x {a \Y b {a} ). 

Remark. The same result x — x^ a ' b ^ = x^ b ' a ^ is obtained exchanging the order of transformations, 
exchanging also the indices a and b in formulas H4.9fH4.13fl . 

Remark. If we denote by x^ 1 ' 2 ^ the quadrilateral lattice obtained by superposition of two (scalar) 
fundamental transforms from x to x^ and x^ 2 \ then the points x, x^\ x^ and x^ 1 ' 2 ^ are 
coplanar again, i.e., the fundamental transformations reproduce the planarity constraint responsible 
for integrability of the quadrilateral lattice. 

4.2. The BQL (Moutard) reduction of the fundamental transformation. In this section 
we describe restrictions on the data of the fundamental transformation in order to preserve the 
reduction from QL to BQL. As usuall (see, for example |211 1151 \'62\ ) the reduction of the funda- 
mental transformation for the special quadrilateral lattices mimics the geometric properties of the 
lattice. Because the basic geometric property of the (scalar) fundamental transformation can be in- 
terpreted as construction of a "new level" of the quadrilateral lattice, then it is natural to define the 
reduced transformation in a similar spirit. Our definition of the BQL reducion of the fundamental 
transformation is therefore based on the following observation. 

Lemma 17. Given B- quadrilateral lattice x : 1 N — > P M and its fundamental transform x con- 
structed under additional assumption that for any point x of the lattice and any pair i,j of different 
directions, the four points x, xi^-s , x/^ and xij\ are coplanar. Then the lattice x : 7L N — > P M is 
B- quadrilateral lattice as well. 

Proof. The result is equivalent to the 4-dimensional consistency of the BQL lattice. Indeed, in 
Proposition let the forth direction be identified withe the transformation direction (the three first 
directions are the lattice directions i, j, k). Then the implication x^ G (£124, £134, £234} is rewritten 
in the form x G (x^ , x {lk) , x (jk) ) . □ 



Definition 2. The fundamental transform x of a B-quadrilateral lattice x : Z w — » P M constructed 
under additional assumption that for any point x of the lattice and any pair i, j of different directions, 
the four points x, Xuf\, xu\ and xq) are coplanar is called the BQL reduction of the fundamental 
transformation of x. 

On the algebraic level, the Darboux-type transformation of the solutions of the linear problem 
12.9fl was introduced and studied by Nimmo and Schief in |38| as discretization of the Moutard 
transformation. We will derive their results from the general theory of transformations of the 
quadrilateral lattice. 

Lemma 18. Given a scalar solution Yi of the linear problem (|2.18(l with the rotation coefficients 
restricted by the BQL reduction (I2.22(l - (|2.23(1 . denote by 6 = Q(Y,h) the corresponding potential, 
where the Lame coefficients are given by equation l|2.21(l . Then the functions 

(4.18) Y* = m > —{6 + 6 

T(i) 

are solutions of the adjoint linear problem H2.16|l in the BQL reduction, and the function 6 2 can be 
taken as the corresponding potential Cl(Y,Y*) 

(4.19) 6 2 = Q(Y,Y*). 

Proof. By direct calculation one verifies that the functions defined in (|4.18|l satisfy the reduced 
system (I2.16|) . Similarly one checks validity of equation i|4.19ll . □ 

Remark. Notice that equation Ij4.19|) implies, under assumptions of Lemma IT8l the form H4.18(l of 
the solution of the adjoint linear problem. 

Proposition 19. Given BQL lattice x with homogeneous representation x in the gauge of the linear 
problem Ij2.9|) . then the transform of x constructed using equation (|4.4|) with the data described in 
Lem,ma \Tfll data satisfies the conditions of the BQL reduction. 

Proof. The fundamental transform of x constructed with such a data reads 

(4.20) x = x- n(X,Y*)/6. 

Equation l]4.2[l . with tl(X , Y*) given above, can be rewritten then in the following form 

Q 

(4.21) X(i) - x = — (x - X(i)) , 

which, together with the linear problem l|2.9l) . implies the linear relation 

6 6 1 

(4.22) — (x (i) - a?) - — (% + (x {ij) -x)=0, i< j, 

C(j) J 

between the homogeneous coordinates of the points x, Xuj) of the lattice and the points xu\ and 
of its fundamental transform, which is the algebraic expression of their coplanarity. □ 

In the approach of [2H1 the Moutard transform x of x was defined in terms of the system Q4.21fl . 
Then x satisfies new Moutard equations (|2.9|l with new potential 

(4.23) pi=f3 e jA, i< j, 

66 {lj) 

and new r-function 

(4.24) f = 6t. 

Another important ingredient of |3H] was the existence of the potential S(6\x) — 6x which satisfies 
the system 

(4.25) A i S(6\x) = 6 ( i ) x-6x {i) . 



We have shown that the algebraic reduction, described in Lemma El of the data of the funda- 
mental transformation can be interpreted as a BQL reduction of the transformation. We close this 
Section by showing that the above algebraic description holds generally. 

Proposition 20. Any BQL -reduction of the fundamental transformation can be algebraically de- 
scribed as in Provosition \lfA 

Proof. We will follow the reasoning of 0] used to the same linear problem ((2.9(1 but in different 
geometric context. Because the BQL-reduced fundamental transformation can be considered as 
construction of the new level of the B-quadrilateral lattice, its algebraic representation should be (in 
appropriate gauge) in the form of the BQL linear problem 12.9(1 

(4.26) x (i) -x = / 0i (x-x (i) ), 

where we can label the transformation direction by the index 0. The compatibility of system (14.26(1 
gives the following equations (compare with ijUnj) 

(4.27) f^f° 3 = f^f\ i-/ ( °/ ) (/ 0i +n = -/(°/ ) r. 

First of them implies the existence of a potential such that 

f 0i = J_ 

thus equations 1(4.21(1 . The second equation rewritten in terms of the potential implies that 8 satisfies 
linear problem l|2.9|l . i.e. 9 = fl(Y,h), where h given by (|2.21l) and Yi is a solution of the linear 
problem l|2.18|l with the rotation coefficients restricted by the BQL reduction (|2.22l) - (|2.23|) . By 
equation (|4.18() we define the corresponding solution of the adjoint linear problem. Finally, direct 
calculation with the help of the Moutard transformation formulas ((4.21(1 show that the potential 

n(X,Y*) = 6(x-x), 

(compare with equation 1(4.20(0 does satisfy equation 1)4.2(1 . thus f2(Y, Y*) is of the form given in 
equation 1(4.19(1 . □ 

4.3. The BQL reduction of the vectorial fundamental transformation. In this Section we 
propose the restrictions on the data of the vectorial fundamental transformation, which are compat- 
ible with the BQL reduction. 

Proposition 21. Given solution Yi : Z N — > V of the linear problem (|2.18|l corresponding to the BQL 
linear problem ((2.9(1 satisfied by the homogeneous coordinates x of the BQL lattice x : 2j N — » P M . 
Denote by & = fi(Y, h) the corresponding potential, which is also new vectorial solution of the BQL 
linear problem 1(2. 9J) . 

1 ) Then 

(4.28) Y* = (-ljEw^'lfe' + eU 

T(i) 

provides a vectorial solution of the adjoint linear problem, and the corresponding potential £l(Y,Y*) 
allows for the following constraint 

(4.29) n(Y,Y*) +n(Y,Y*f = 20® 0*. 

2) The fundamental vectorial transform x of x, given by (14.4(1 with the potentials ft restricted as 
above can be considered as the superposition of dimV (scalar) discrete BQL reduced fundamental 
transforms. 

Proof. The point 1) can be checked by direct calculation. To prove the point 2) notice that when 
dim V = 1 we obtain the BQL reduction of the fundamental transformation in the setting of Propo- 
sition E] For dimV > 1 the statement follows from the standard reasoning applied to superposition 
of two reduced vectorial fundamental transformations (compare with |2"THl"5] V 



Assume the splitting V = V a © of the vectorial space V, and the induced splitting of the 
basic data Y i of the transformation. Then we have also (in shorthand notation, compare equations 



(4.30) 

and the constraint 
(4.31) 



& = 

reads 

n a b 



& a 
& b 



n(Y,Y*) = 



n b n b 



nr n 



e° ® © at © a ® 



0"(g)0 



0"(g)0 



By straightforward algebra, using equations l|4.31|l . one checks that the transformed potentials (com- 
pare equations 14.13J1 ') 



(4.32) n b b {a} = ft b b - 

(4.33) b{a} = b 
satisfy the BQL constraint 14.29(1 as well, i.e., 



5 L bi 



■ n^n-]- 1 ©- 



(4.34) 

which concludes the proof. 



b{a} 



b{a}t =2&b{a} Q b{a}t 



□ 



Remark. Because the BQL-reduced fundamental transformation can be cosidered as construction 
of new levels of the B-quadrilateral lattice, then if we denote by x^ 1,2 ^ the B-quadrilateral lattice 
obtained by superposition of two (scalar) such transforms from x to x^ and x^ 2 ' , then for each 



direction i of the lattice the points x, x 



{1} ^{2} 



and x^ 1 ' 2 } are coplanar as well as the points x 



(0' 

x x ^ J , x^ and £^ 2 \ Similarly, if we consider superpositions of three (scalar) transforms of the B- 
quadrilateral lattice x then the points x, x''- 1,2 '', x^ 1 ' 3 ^ and x^ 2 ' 3 ^ are coplanar as well as the points 

£{2} } £{3} and £{1,2,3} 

4.4. The Pfaffian form of the transformation. Finally, we are going to show that the Pfaffian 
formulas of the vectorial discrete Moutard transformation obtained in [38] can be derived from the 
corresponding formulas of the fundamental transformation subjected to the BQL reduction. 

Denote by S(0|0) the antisymmetric part of £l(Y,Y*) then equations (|4"TJl . 62B) and 
imply 



(4.35) 



A,S"(0|0) = © (i) <g> 0* - ® 








-0* 
S{®\®) 



Lemma 22. For Cl(Y,Y*) and 5(0|0) as above we have 

(4.36) det n(Y", Y*) = det 5(0|0) 

Proof. Notice that the jth column ttj of tt(Y, Y*) is of the form 

(4.37) Slj = i e + S j) 

where 9^ is the jth component of 0, and Sj is the jth column of 5(0 10). Then the basic properties 
of determinants imply that 

(4.38) 



dim^ 



dctn(Y,Y*) =dctS(0|0)+ > 9 J S(j), 



where by S(j) we denote the matrix 5(0|0) with jth column replaced by 0. The second summand 
in (|4.38|) is the Laplace expansion of that in l|4.36l) . □ 

The standard properties of determinants of antisymmetric matrices (see Appendix [5J imply the 
following result derived in [3HJ directly on the level of vectorial Moutard transformation. 



Corollary 23. The transformation formula ()4.5|l of the QL T-function and the relation (|2.25[1 
between both T-functions imply the following transformation formula for the BQL T-function 

'rPf 5(0|S), 

-0*, \ 
5(010) ) ' 



(4.39) 



Pf 



dimV even 
dimV odd. 



Remark. Notice jHH] that 

-0*, 
5(01©) 

which allows to define 



(4.40) 



S(0|0), where 



1 





(4.41) 

and gives 
(4.42) 



V{&) 



Pf 5(0|0), 
Pf 5(010), 

f = tV(&) 



dim V even, 
dim V odd, 



Finally, we will connect the formula of the vectorial fundamental transformation l|4.4|l in the BQL 
reduction with the Pffafian form of the vectorial Moutard transformation 38 . 

Corollary 24. The homogeneous coordinates (in the gauge of the linear problem ()2.9(l ) of the BQL 
lattice x obtained from the BQL lattice x via vectorial transform with the solution of the linear 
problem (|2.9|l are given by 

Ai P(®,x*) 



(4.43) 
where 
(4.44) 



V(@) 



T{®,x l ) = V 





x i 



Proof. We will work using assumptions and notation of Proposition Let us define 

(4.45) S{x\9)=to{X,Y*)-x®&, 
then, due to equations l|2.17jl . (|4.2|) and (|4.28|l we have 

(4.46) AjS^Jcle) = cc (i) (8 0* - x ® 0^. 

By the Cramer rule and equation l|4.37(l . formula 1|4.4JI in the considered reduction case can be 
brought to the form 

/ detS'(l) 

x = x - (x «) 0* + S(x\&)) 



(4.47) 
moreover we have 
(4.48) 



1 



dctn(Y,Y* 



\ det5(dimV) 



/ detS(l) 



l 



-0* 
5(010) 



\ detS(dimV) 

Our further analysis splits in the cases of dimV being even or odd. In the first case the right 
hand side of equation 14.48|) vanishes giving 

PfS[l] 



(4.49) 



x = x — 



S{x\&) 
PfS(©|0) 



Pf5[di 



(4.50) 



1 



where we used equation (|4.36|) and the Pfafhan analogue (|B . 1 1|> of the Cramer rule for for solutions 
of the equation S(®\&)y = 0. Then the expansion rule for Pfaffians 1|B.8(1 implies that the ith 
coordinates of the B-quadrilateral lattice x and its transform x can be put in the form 

x l S(at\8) 

-x l -0* 

-s{x i \®y © s(0|©) 

For dimV odd, by equations l|4.36|l and 14.48(1 . equation (|4.47() reduces to 

! / detS(l) 



PfS(0|0) 



-Pf 



g(0_y 
V(&) 



(4.51) 



x = — 



-0* 
S(&\®) 



S(x\@) 



\ detS(dimV) 



Expanding det S(j) with respect to its jth column, and using Pfafhan expressions (|B.13J) for the 
minors of ^(©j©) we obtain 



Pf 



(4.52) x 
which concludes the proof. 



S(x l \&) 
-5(^1©)* S(0|0) 



Pf 



-0* 
S(®\&) 



V{&) ' 



□ 



5. Conclusion and remarks 

We presented new geometric interpretation of the discrete BKP equation within the theory of 
quadrilateral lattices. This new integrable lattice should be considered, together with the symmet- 
ric lattice |18| and the quadrilateral lattices subject to quadratic constraints |15j . as one of basic 
reductions of the quadrilateral lattice. In the forthcomong paper |lfi| we show, for example, that 
the discrete isothermic surfaces are lattices subjected simultaneously to the BQL and quadratic 
(in this case the quadric is the Mobius sphere) reductions. 

As in the case of the Hirota (the discrete KP) equation, also the Miwa (the discrete BKP) equation 
can be considered in the finite fields (or the finite geometry) setting. In paricular, the main algebro- 
geometric way of reasoning (see 20, 5] for the former discrete KP case) leading to the Prym varieties 
should be also transferable for fields of characteristic different from two (see [HZ] for general theory 
of Prym varieties) . 

The explicit Prym-theta functional formalae (it is enough to consider the case N = 2) for the wave 
function and potentials of the discrete Moutard equation can be also used to provide characterization 
of the Prym varieties among all principally polarized abelian varieties (the Prym-Schottky problem) 
in the spirit of 
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Appendix A. An alternative proof of the existence of the BQL gauge 

The planarity condition of elementary quadrilaterals of QL can be expressed in terms of generic 
homogoneous representation as the following system of discrete Laplace equations |17| 

(A.l) X(ij) = a tj X(i) - a ri £C(j) + c lj x, 1 < i < j < N 



whose compatibility are equations 

( A 9\ n i] r lk - n jl r^ k — n ik — n ki r^ k — n jk - n kj r ik 

(h -i\ J3 „»* — n ij — J k i J k n V _ n kj n ik 

[A.6) a (k) a - a {j) a - + a (i) a a {i) a 

( k A\ ni k — n^ k — —r ik 4- n ik 4- n ki n^ k 

lA ^ n¥ n kl — n kj n kt — r i] 4- n ji n kl - n i3 n U 

where l<i<j<k<N. Because 

(A.6) x A !E(y) A a; (ife) A x (jk ) = x A A ccq-) A x (k ^(a lj a jk a kl - a Jl a kj a lk ), 

then the BQL reduction condition is equivalent to 

(A.7) a zj a jk a kt - a Jl a k] a tk = 0. 

We will show that equation (jA. T|> implies existence the gauge function p : T, N — > R such that 

(A.8) a ij p (i) = a ji p u) , a lk p {i) = a kl p [k) , a lk p {j) = a k] p [k) , 

(A.9) p (tj) = c l V, p (lfc) = c ik p, p (]k) = c jk p. 

Then, after rescaling x — > x/p, the new homogeneous coordinates satisfy the system H2.9[) . 

Let us consider equations (|A.8I) - (|A.9(I as a difference system, which allows to calculate from p and 
(say) values of the gauge function in remainig vertices of the hexahedron. Notice first that the 
condition (|A.7|I and the system l|A.2() - ilA.5|) imply (it can be verified directly, but actually it follows 
from Corollary |5l that 

/ a 1 n\ *J ik ji ik ik ij ki ik jk ij kj ik 

(A.10) a { J k) c =a J [k) d> , a (]) c J = a {]) c> , a J (i) c = afic . 

The condition (|A.7(l assures self-consistency of equations (|A.8|) . Then equations IjA.lOfl imply 
consistency of equations (|A.9fl with the following consequence of equations 



a \k)Pm = a ll)P(jk), a l (j)P(i 3 ) = a u)P(jk), a^Piij) = a^Pftk)- 



Finally, the self-consistency of equations ljA.91) in finding pujk) follows from equations (| A. 8|> and the 
system (|A.3ll - (|A.5|l . 

Remark. Notice that the system l|A.2j) - (|A.5|) of nonlinear equations can be considered as a system 
of eight linear equations allowing for transition 

(n^ n ik n ki n^ k n k 3 r ik r^ k \ I n ij n ji n ik n ki n jk n kj r ij r ik r^ k \- 

(a ,a ,a ,a ,a ,a ,c ,c ,c ) — >• (a (fe) , a (fe) , a^, a^, a (i) , a (i) , c (fe) , c^, C (i) j, 

the difference in the number of unknowns and equations reflects the homogeneous nature of the 
linear system QA.1[1 . 

Appendix B. Pfaffians 

We recall basic properties of Pfaffians [SHE], which we use in Section POI Let A = (<Xij)i<i,j<2r 
be a skew symmetric matrix (i.e., a,ji = —dij) of the even order 2r. Consider the form 

(B.l) uj = ajjej A e 3 , 

i<j 

then the Pfaffian Pi (A) of A is defined by 

(B.2) u Ar = (r!) Pi(A) e x A ■ ■ ■ A e 2r . 

For each permutation ir of {1, . . . , 2r}, put A* = (fl,^,^), then 
(B.3) Pi{A") =sgn7rPf(A). 



Notice the analogy with the determinant det£? of an arbitrary square matrix B — (6y)i<j,j< n 
expressed in terms of the forms 

(B.4) uJi = J2 b ii e j 

3 

as follows 

(B.5) u>i A ■ • • A uj n — det(-B) e\ A • • • A e n . 

It turns out that the determinant of any skew symmetric matrix of an even order equals to the 
square of its Pfafhan 

(B.6) det(A) = (Pf(A)) 2 . 

For any two subsets I,JC {1, . . . , n} denote by B(I, J) the sub-matrix of B obtained by removing 
all the zthE I rows and all the jthe J columns of B. Then in analogy to the Laplace expansion of 
determinants 

n 

(B.7) % det(B) = hj(-l) k+i det(B({fc}, {*})). 

fc=i 

we have the following expansion formula for Pfafhans 

2r 

(B.8) SijPfiA) =^a fc ,(-l) fc+i - 1 Pf(A({fc,z},{fc, ? })). 

k=l 

Both formulas imply 

(B.9) det(A({i}, {j})) = -Pf(A) Pf(A({i, j}, 

which leads to the following Pfaffian-Cramer rule for solutions of the linear system 

(B.10) Ay = b 

with non-degenerate skew symmetrix matrix od the even order: 



(B.H) V - p{(A) , 

where by A\j] is denoted the matrix A whose jth column is replaced by 

(B.12) b'=(b\...,V-\0,-V+\...,-b 2r )\ 

and whose jth row is replaced by — b' . 

When the order of the skew symmetric matrix A is odd we have det(A) = 0, but the following 
formula holds 

(B.13) detL4(«, {j})) = P{(A({i}, {<})) Pi(A({j}, {j})). 
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